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Abstract 

We present the results of Monte Carlo simulations of the magnetic properties of individual spherical nanoparticles 
with the aim to explain the role played by surface anisotropy on their low temperature magnetization processes. 
Phase diagrams for the equilibrium configurations have been obtained, showing a change from quasi-uniform 
magnetization state to a state with hedgehog-like structures at the surface as ks increases. Through the simulated 
hysteresis loops and the analysis of spin configurations along them, we have identified a change in the magnetization 
reversal mechanism from quasi- uniform rotation at low ks values, to a non-uniform switching process at high ks- 
Results for the dependence of the coercive field and remanence on ks and particle size are also reported. 
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1. Introduction 

When reducing the size of magnetic particles to 
the nanoscale, surface and finite-size effects play a 
dominant role on establishing their peculiar low tem- 
perature magnetic properties [1]. Reduction of the 
magnetic ordering temperature and spontaneous mag- 
netization as particle size is reduced [2,3], increased co- 
ercivities and high closure fields in the hysteresis loops 
[4] are often attributed to the frustration associated 
to broken bonds and disorder present at the particle 
surface. However, in real samples, these effects are 
usually masked by the presence of particle size distri- 
butions and by the magnetic interactions between the 
particles (often mentioned as a possible source of their 
spin-glass-like behavior), making it very difficult to 
discern what is the real contribution of surface effects 
alone on the magnetic properties. Moreover, current 
experimental techniques do not allow to probe sur- 
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face microscopic structure in geometries such as those 
of a nanoparticle, a fact that forces to assume phe- 
nomenological models for surface anisotropy. Monte 
Carlo (MC) simulations constitute a benchmark to 
probe different models since, by modeling magnetic 
ions as classical [5,6,7,8] or Ising spins [9,10,11], the in- 
teratomic interactions, the local magnetic anisotropy 
directions and the values of anisotropy constants can 
be easily varied while taking into account the exact 
geometry of the underlying lattice of magnetic ions. 
Here, we will present results of MC simulations aimed 
to understand the influence of surface anisotropy on 
the magnetic properties and hysteresis of a single 
nanoparticle, neglecting the effects of interactions with 
other particles. 



2. Model and Results 

Since most of the peculiar low temperature magnetic 
properties of nanoparticles are more clearly observed 
in ferromagnetic oxides, we will consider here a parti- 
cle with lattice geometry of maghemite (7-Fe2 0a), a 
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ferrimagnetic compound in which the magnetic ions sit 
on two sublattices (O, T) with different coordination. 
We have considered particles with spherical shapes of 
diameters D — 3 — 6a (a is the lattice constant) .The N 
magnetic ions will be modeled by Heisenberg classical 
spins Si of unit magnitude with the following interac- 
tion Hamiltonian in temperature units 

JV 

H/k B = -J2 J v s * s J-J2 h Si 

N c N s 

-J2[k c (S?f] +J2 [ksiSi-m) 2 ] . (1) 

i=l i=l 

This is an extension to Heisenberg spins of our previ- 
ous model for Ising spins in which the first term ac- 
counts for the n.n. exchange interaction, with Jij the 
(intra)inter-sublattice (negative) exchange constants 
(see [11] for numerical details) and the second repre- 
sents the interaction with the external field h (this can 
be translated to physical units through the relation 
H = kshj 'fi, with fj, the magnetic moment of Fe ions). 
The last two terms have been introduced to differenti- 
ate the magnetic anisotropy at the Nc core spins, kc = 
1 (uniaxial along the z-axis) , from that at the Ns sur- 
face spins (those in the outer unit cell), ks, which will 
be considered to point along the local radial direction 
hi due to breaking of crystalline symmetry at parti- 
cle boundaries. Our goal is to understand the change 
in the magnetic properties as ks is increased from kc 
for different particle diameters D. Some of the details 
of the MC method as implemented in our simulations 
can be found in previous works [12,13]. 




Fig. 1. (Color online) Phase diagrams showing the con- 
tributions of (a) the surface and (b) the core of particles 
with diameters D = 3a, 5a to the projection along the lo- 
cal easy-axes M n as a function of the surface anisotropy 
constant ks/kc- The Inset in (a) shows the total magneti- 
zation along the z-axis M z . Below: snapshots of the corre- 
sponding T = configurations of a spherical particle with 
D = 5a obtained after cooling from a disordered state at 
high T are displayed for representative values of ks in the 
two regimes. 

In order to characterize the influence of the surface 
anisotropy on the magnetic order at low temperatures, 



we have first performed simulations to find the equilib- 
rium configurations at low T. For this purpose, we start 
from a disordered configuration of spins at a high tem- 
perature (T = 50K) and follow a progressive cooling 
down to T = in constant temperature steps S — —IK 
in absence of magnetic field. In order to better charac- 
terize the magnetic ordering of the states achieved at 
T = 0, we have also computed the quantity M^' c = 

5Zi=fi C \&* " "»l> which indicates the departure of spins 
from alignment along the local easy-axis (M n — 1 in 
this case). In Fig. 1, we present the values of Mf' c and 
M^' c of the configurations corresponding to ks/kc 
ranging from 1 to 100 for particles with diameters D — 
3a, 5a. The curves in this figure, indicate ferrimagnetic 
order along the z-axis when M% ~ 1 and M„ ~ 0.5 and 
radial alignment of spins when Mn ~ 0.5 and M« ~ 1. 
As can be clearly seen, at a value ks = k s , there is a 
change in the magnetic order of the configurations as 
marked by the sharp decrease of M% and \M Z \ to zero 
and the increase of M® to 1. For ks < k s the particle 
has a ferrimagnetic ordered core along the z-axis (see 
left upper snapshot in Fig. 1) with surface spins that 
progressively deviate towards the radial direction as ks 
increases. This gives rise to the formation of vortex- 
like structures that extend to the particle core, as can 
be seen in the right upper snapshot in Fig. 1. For ks > 
k s , the net magnetization of the particle tends towards 
zero. This is due to the fact that the department of sur- 
face spins towards the radial direction is transmitted 
to the core via the antiferromagnetic exchange interac- 
tions, giving rise to hedgehog-like structures (see ower 
snapshots in Fig. 1). 
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Fig. 2. (Color online) Hysteresis loops for a particle with 
D = 5a for different values of ks as indicated. The differ- 
ent panels stand for (a) the total magnetization, (b) the 
contribution of the core to M n , and (c) the contribution of 
the surface to M n . 

We have also simulated the hysteresis loops of par- 
ticles with diameters D — 3a, 5a for different values 
of ks- We initially apply a strong enough magnetic 
field along the z-axis strong that is subsequently var- 
ied in steps Sh — 1 K during which the magnetization 
is averaged over 1000 MC steps at each field value. In 
Fig. 2, we present hysteresis loops for a particle with 
D — 5a for representative values of ks- Let us first no- 



tice that, as surface anisotropy increases, the hystere- 
sis loops change from quasi-squared shape, indicating 
that the dynamics of the particle may be described as 
the quasi-uniform rotation of the magnetization due 
only to the non-compensated spins, to more elongated 
shape with higher closure fields as observed in oxide 
particle systems [14] . This is in agreement with the for- 
mation of disorder states at the particle surface due to 
increased anisotropy and frustration due to finite-size 
effects. This change in the reversal mechanism can be 
more clearly analyzed by comparing panels (b) and (c) 
in Fig. 2, where, for values of ks > k s , we see that sur- 
face spins reverse by forming radial structures (M® ~ 
1 along the loop) that induce the subsequent reversal 
of the core [15] . Further insight into the reversal mech- 
anism can be gained from the dependence of the coer- 
cive field he and remanent magnetization Mn em on ks 
presented in Fig. 3. The increase of he and the decrease 
in M_R em observed in Fig. 3, are clear indications of a 
change in the magnetization reversal mechanism from 
quasi-coherent to a surface dominated process [12] . The 
snapshots of the remanent configurations shown in Fig. 
3 demonstrate the formation of a shell with disordered 
spins pointing towards the radial direction which in- 
creases in thickness as ks is increased. This last feature 
also explains the obtained dependence of he and M_R em 
on ks- Moreover, he increases with decreasing size as 
observed by Tronc and coworkers in experiments with 
maghemite particle dispersions [16,17]. 
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Fig. 3. (Color online) Dependence of (a) the coercive field 
he and (b) remanent magnetization Mn em on % as ob- 
tained from low T hysteresis loops for spherical particles 
with diameters D = 3a, 5a. Lower panels: remanent states 
along the hysteresis loop of the D = 5a particle at represen- 
tative values of ks in the two reversal regimes, snapshots 
display spin configurations at the equatorial plane perpen- 
dicular to the z-axis. 
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